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ABSTRACT 

The SoHO spacecraft made astrophysicists achieve a major breakthrough in 
the knowledge of the Sun. In helioseismology, both GOLF and MDI experiments 
aboard SoHO greatly improve the accuracy of seismic data. More specifically, 
the detection of an enhanced number of low degree low order modes improves 
the accuracy on the sound speed and density profiles in the solar core. After a 
description of these profiles, we build solar seismic models. Different models are 
considered and enable us to derive precise emitted neutrino fluxes. These ones 
are validated by the seismic data and are in agreement with the recent detected 
neutrinos, assuming 3 neutrino flavors. The seismic models are also used to put 
limits on large scale magnetic fields in the solar interior. This analysis puts 
some upper bounds of ~ 3 x 10 7 G in the radiative zone. Such a field could 
slightly improve the emitted neutrino flux, which remains in agreement with 
the Sudbury Neutrino Observatory result of 2001. From the models we deduce 
gravity mode predictions, and the electron and neutron radial densities that are 
useful to calculate the neutrino oscillations. We also begin to discuss how the 
external magnetic field may influence such quantities. 

Subject headings: instrumentation: GOLF, MDI — physics: neutrinos, magnetic 
fields — spacecraft: SoHO — Sun: interior, oscillations 



1. INTRODUCTION 

Since the detection of oscillations at the surface of the Sun in the early sixties (Leighton, 
Noyes & Simon 1962), and their later interpretation as trapped acoustic waves, helioseismol- 
ogy turned out to be a very powerful tool to probe the solar interior. Compared to ground 
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networks, the Solar and Heliospheric Observatory spacecraft (SoHO) allows very long and 
continuous measurements leading to very low amplitude detections (down to 3 mms -1 after 4 
years): we detect the low frequency part of the oscillation spectrum, not polluted by the tur- 
bulent surface. We also gained confidence in the frequency measurements from comparisons 
between both the Global Oscillation at Low Frequency (GOLF, see Gabriel et al. 1995) and 
Michelson Doppler Interferometer (SOI/MDI, see Scherrer et al. 1995) instruments aboard 
SoHO: these comparisons pointed out the consistency between the two data sets for modes 
between 1.4 and 3.7 mHz (Bertello et al. 2000a). The SoHO data are of high quality, since 
roughly speaking the frequencies measured with GOLF are affected by an error one order 
of magnitude smaller than the one of 1995. Moreover, the data quality was enforced by 
progresses in the spectral analysis methods. For instance, Bertello et al. (2000a) use the 
new RLSCSA method, while Garcia et al. (2001) use probability computations and different 
spectral methods (periodogram, random-lag average cross-spectrum, homomorphic decon- 
volution...). The improvement of seismic data, the addition of new low degree low order 
p modes (Bertello et al. 2000b; Garcia et al. 2001), and the comparison between indepen- 
dent instruments make us confident in the results of the inversions based on the GOLF and 
SOI/MDI data. 

Therefore, we utilize these inversions to build solar seismic models. Usually, the term 
"seismic models" refers to models that have been directly derived from seismic data (e.g. 
Kosovichev & Fedorova 1991; Dziembowski et al. 1995; Basu & Thompson 1996; Shibahashi 
& Takata 1996; Antia & Chitre 1998). These models use the primary inversion of seismic 
data that returns the sound speed c s (r) and the density p(r) profiles inside the Sun. By 
assuming hydrostatic equilibrium, it is straightforward to derive the pressure profile P(r). 
With additional conditions on the input physics, the temperature T{r) and helium abundance 
Y{r) profiles can be determined. The knowledge of these different quantities as a function 
of the fractional radius defines a "seismic model" . 

In our case this term refers to something different: we stay in the classical framework of stellar 
evolution and we use a ID stellar evolution code to compute a model in agreement with the 
seismic observations. This agreement is obtained by varying a few physical parameters inside 
their error bars. This is possible now, thanks to 10 years of improvements in solar modelling 
with the introduction of updated physics, microscopic diffusion, & turbulence at the base of 
the convection zone (e.g. Turck-Chieze et al. 1988; Turck-Chieze & Lopes 1993; Dzitko et 
al. 1995; Brun, Turck-Chieze, & Morel 1998; Brun, Turck-Chieze, & Zahn 1999). 

The seismic models detailed in this paper are used to derive the neutrino fluxes emitted 
by the Sun. The main interest is that these fluxes include observational data, meaning that 
helioseismology can help in properly determining the neutrino production: a key point for 
the neutrino puzzle. 
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We have presented this first result in Turck-Chieze et al. (2001b) showing a remarkable 
agreement with the detections of the Sudbury Neutrino Observatory (hereafter SNO) plus 
SuperKamiokande (SK) experiments. In this paper, we detail and generalize our approach. 
In section 2 we introduce the data we use and the major features of the evolution code that 
permits us to compute models. In section 3 we discuss the way the sound speed and density 
inversions are carried out. In section 4 we derive the seismic models and discuss about their 
uniqueness. In section 5 we compute the seismic oscillation frequencies. In section 6 we 
compare the predicted neutrino fluxes to the ones detected by the terrestrial experiments 
like SK or SNO, and derive some quantities useful to the determination of the neutrino 
oscillation parameters. In section 7 we focus on the magnetic field and the way we can put 
upper bounds on it. Finally we conclude in section 8. 



2. THE HELIOSEISMIC OBSERVATIONS AND THE PHYSICAL INPUTS 

OF THE SOLAR MODELS 

To derive the sound speed (c s ) and density (p) of the Sun, we use the data from both 
the GOLF and MDI instruments and the inversion procedure described in the next section. 
For the p modes with i < 2, we use Bertello et al. (2000a, 2000b). These frequencies are 
extracted from GOLF data. For the p modes with £ > 3, we use Rhodes et al. (1997) who 
utilized data from the SOI/MDI instrument. The quality of these modes has been discussed 
in Turck-Chieze et al. (2001b) showing the results of some inversions with different values of 
the radial order n; in section 5 we show, by comparing observed and calculated frequencies, 
how the SoHO observations have improved the quality of the seismic data. 

The solar models are computed with the CESAM code (Code d'Evolution Stellaire 
Adaptatif et Modulaire, see Morel 1997). This is a ID quasi-static stellar evolution code 
that solves the stellar structure equations by a spline collocation method. We always start 
the evolution from the pre- main sequence (PMS). The basic physical characteristics of the 
models are: 

• the nuclear reaction rates from Adelberger et al. (1998), with Mitler intermediate 
screening (Mitler 1977). For the 7 Be(p,7) 8 B reaction, we use the S(0)n value derived 
by Hammache et al. (1998). For the 7 Li(p, 4 He) 4 He astrophysical factor we use Engstler 
et al. (1992). In order to properly compute the lithium burning on PMS, we adjust 
the time step and the rotation law, according to Piau & Turck-Chieze (2001); 

• the opacities are derived from the OPAL95 opacity tables (Iglesias & Rogers 1996) 
for temperatures larger than 5600fC For lower temperatures, we use the Alexander 
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opacities (Alexander & Ferguson 1994); 

• the equation of state (EOS) is OPAL (Rogers & Iglesias 1996); 

• the microscopic diffusion is taken into account with the prescription of Michaud & 
Proffitt (1993); 

• the turbulent mixing at the base of the convection zone (hereafter BCZ) is treated, 
following Brun et al. (1999). 

The references listed emphasize the improvement in the solar interior physics that oc- 
cured in the last decade, in parallel to the improvement of the seismic data. This allowed us 
to reject some extra physical processes like large screening and mixing in the central region 
(e.g. Turck-Chieze et al. 2001a). 

3. SEISMIC DATA AND THE INVERSION OF THE SOLAR SOUND 

SPEED AND DENSITY 

The inversion results for the sound-speed and density profiles are obtained by using the 
Optimally Localized Averaging (OLA) method (see Kosovichev 1999 — hereafter K1999— 
for mathematical details). In this method the differencies between the observed and model 
frequencies are expressed as a sum of two linear integrals for the corresponding relative 
differences in the sound speed and density (see eq. [52] of K1999). The sensitivity kernels 
in these integrals are calculated by using a variational principle for adiabatic non-radial 
oscillations that allows us to neglect the variations of the eigenf unctions to the first order 
of approximation. The non adiabatic effects, and the uncertainties in the physics of the 
near-surface layers, are taken into account by adding a term to the two integrals related to 
the frequency differences. This term is a smooth function of frequency, scaled with the mode 
inertia (see eq. [72] of K1999). 

The estimates of the localized averages for the sound-speed and density corrections to 
the solar models are obtained by considering linear combinations of the integral relations for 
the frequency differences. We proceed such that the corresponding linear combinations of the 
sensitivity kernels form narrow localized, Gaussian-type, kernels at various target positions 
along the solar radius for one of the variables (sound speed or density) and are negligible for 
the other variable. The inversion procedure includes additional constraints to eliminate the 
surface term approximated by a linear combination of Legendre polynomials of degree less 
than 5, and also to minimize the errors of the sound-speed and density corrections. The later 
constraint is based on observational error estimates of the mode frequencies, and includes 
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a regularization parameter that controls the trade-off between the spatial resolution of the 
inversions (measured as a "spread" of the averaging kernels) and the error magnification. 
The regularization parameter is chosen to provide a sufficiently smooth radial dependence 
of the sound-speed and density corrections. 

The results are presented in our figures in the form of the horizontal and vertical bars 
centered at the center of gravity of the localized averaging kernels. The size of the horizontal 
bars corresponds to a characteristic width ( "spread" ) of the averaging kernels and provides 
an estimate of the spatial resolution, and the vertical bars correspond to la formal error of 
the sound-speed and density corrections (see Table 1 of Turck-Chieze et al. 2001b). The 
standard solar model of Brun et al. (1999) was used as a reference. The inversion procedure 
was tested by using various other solar models as Sun's proxy and adding random Gaussian 
noise to calculated frequencies of these models to simulate the observational errors. 

Figure 1 shows the inversions obtained for the sound speed and the density. This figure 
highlights how the low order modes impact deeply on the density profile and how they 
slightly improve the sound speed in the solar core. In the following, we shall call seismic 
model a model that reproduces as well as possible the observed sound speed profile. We will 
not take the density as a reference yet, since we first wait for confirmation of its profile by 
new detections at low frequencies, due to its extreme sensitivity to these modes. The density 
appears more sensitive to the detailed physics than the sound speed. It will be used in the 
future to progress on the dynamics of the Sun. 

4. BUILDING THE SEISMIC MODELS 

4.1. Starting Point 

In this study, our goal is to cancel the discrepancy on c s between the solar models and 
the real Sun. To this end, we slightly modify a few physical parameters used in CESAM. Of 
course, we derive seismic models that are not unique: we could obtain similar results on the 
sound speed by adjusting different parameters. However, we justify our choices in section 
4.3. We adjust the physical quantities the sound speed is sensitive enough to, and modify 
them within their error bars. It is decided to change as little parameters as possible. 

We start from a specific solar model: the non-standard Btz model of Brun et al. (1999). 
This model was designed to reduce the discrepancy on c s between the standard models and 
the Sun at the BCZ, by taking into account the horizontal motion produced by the sudden 
disappearance of the differential rotation profile in the tachocline. Doing so, they get for the 
first time a correct 7 Li abundance at the solar surface (Ar Li = 1.16 ± 0.10 dex according to 
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Grevesse & Noels 1993 [hereafter G&N93]): the only abundance that was not predictable 
before. It is a real improvement since the standard solar models do not consider any turbulent 
mixing at the BCZ. It is well known that the lithium is the only indicator of the internal 
structure for numerous stars and its abundance is very difficult to predict in classical stellar 
evolution. However, the adopted tachocline prescription is purely hydrodynamic and does 
not account for any magnetic field, eventhough the magnetic dynamo process is thought to 
occur in the thin tachocline. 

The main impacts of this prescription are: to reduce the influence of the microscopic 
diffusion — the diffusion of heavy elements toward the center of the Sun is slowed down — 
and to burn 7 Li on the main sequence. Three parameters define the tachocline: 

• its current width, taken as d — O.O5i? ; 

• the Brunt- Vaisala frequency at the BCZ, N = 25 /zHz; 

• the present rotation rate at the BCZ, fi — 415 nHz. 

Moreover, the turbulent diffusion coefficient in the tachocline is time-dependent, since it is 
related to the rotation rate. In the Btz model the rotation as a function of time follows the 
Skumanich's law (Skumanich 1972). 



4.2. The First Seismic Model: Seismicx 

Starting from Btz, we undertake the construction of a first seismic model (already 
presented in Turck-Chieze et al. 2001b). We focus primarily on the solar core where the 
energy generation and neutrino production occur: to adjust the physical parameters, we 
take advantage of our knowledge of the sensitivity of the model to the physical ingredients 
(through, for instance, Turck-Chieze & Lopes 1993; Dzitko et al. 1995; Turck-Chieze et al. 
1997; Turck-Chieze et al. 2001a). 

We increase the p — p cross section (hereafter S(0) p - P ) by 1%: usually S(0) p - P = 
4.00(1 ± 0.007±omi)10- 22 keV b (Adelberger et al. 1998). As this cross section is only 
known theoretically, we try to constrain it by helioseismology. The change of S(0) p - P has a 
rather large impact on c s . An increase of this cross section induces a decrease of the core 
temperature, since the model is calibrated to obtain the solar luminosity. 

In parallel to this modification, we adjust the initial metallicity (Zq) of the Sun. Com- 
pared to Btz (Zq = 0.01959), we increase it by 3.9%. In the radiative interior, a major 
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contribution to the opacities comes from the heavy elements (see Fig. 2): at O.6i? , about 
85% of the opacities are due to the metals. In the solar core, this part is still 45%, in this 
case all the elements are almost fully ionised, except iron. Actually, the opacities in the 
core are primarily due to inverse brehmstrahlung and electron scattering, while at O.6i? the 
bound-bound processes dominate. By increasing Zq, we change both the Rosseland opacities 
(K ross ) and the mean molecular weight (//), explaining why the sound speed of the seismic 
model decreases in the core and goes up in the radiative zone: in the core, the raise of /i 
dominates the raise in temperature produced by the increase of K ross ] beyond O.3i? , the 
raise in temperature overcomes the increase of /i. A last change is carried out: the OPAL 
EOS is tabulated and depends on Z . Instead of Z = 0.01959 used to compute the EOS 
table for Btz, we use a table computed with Zq = 0.0203. Actually, this change has no 
important impact on the sound speed profile. 

Having reduced the discrepancy on c s in the radiative zone, we undertake to diminish 
it at the BCZ. To this end we modify the parameters defining the tachocline: we reduce its 
width from O.O5i? to O.O25i? , in accordance with the most recent helioseismic results (e.g. 
Elliott et al. (1998) announce d = O.O2i? , Corbard et al. (1999) announce d < O.O5i? ); 
we also increase the present rotation rate f2 from 415 to 430 nHz (Corbard et al. 1999). 
Moreover, to retrieve a good 7 Li abundance at the surface, we increase the Brunt- Vaisala 
frequency N from 25 to 105 /zHz. This frequency undergoes a dramatic change when ap- 
proaching the convection zone. An increase of N reduces the efficiency of the turbulent 
mixing and diminishes the 7 Li depletion. 

Another improvement is the change in the rotation law. Instead of the Skumanich's 
law, which is far less realistic on the PMS, we utilize the law proposed by Bouvier et al. 
(1997), following Piau & Turck-Chieze (2001): the rotation rate of the Sun slightly goes 
down during the early evolution phase, and experiences a rapid acceleration at 10 million 
years when the circumstellar accretion disk separates from the young Sun. Finally, the 
rotation rate decreases following the Skumanich's law because of the magnetic breaking. 
Eventhough the impact of such a rotation law on the solar evolution is rather weak, its use 
is more realistic. To reach a correct 7 Li content at 4.6 Gyr, we suppose that the Sun was a 
slow rotator on the PMS (this is why we choose a disk separation at 10 million years). 

To complete the building of the seismic model, we concentrate on the outer part of the 
Sc 2 /c 2 profile: it is still far from flat, since no appropriate model exists for the upper layers 
with a ID code. However, it is possible to obtain a better agreement with the Sun by cali- 
brating the seismic model at a radius R\ different from the standard one R & = 6.9599 x 10 10 
cm deduced from photometric observations (e.g. Allen 1976). Recent analyses based on f- 
mode frequencies (e.g. Schou et al. 1997; Antia 1998) worked out seismic radii, respectively 



- 8- 



R = 695.78 Mm and R = 695.68 Mm, slightly smaller than the photometric one. This latter 
is greater than the recent optical determination R = 695.51 Mm by Brown & Christensen- 
Dalsgaard (1998). There is neither a firm answer explaining the origin of these discrepancies, 
nor the effect of the solar cycle on the radius. Thus, we calibrate the seismici model with 
Ri = 6.95936 x 10 10 cm, with a consequent ad hoc improvement in the convective zone. 
Nevertheless, the inversion of the density and sound speed was carried out assuming the 
standard value for R Q . 

The Table 1 lists all the main features of the models described in this paper. The seismici 
model is also available on the web site http://apc-p7.org/Neutrino-APC/Sismic-model.html, 
with the detailed values of many parameters for a large number of shells (including the 
electron number density, see sect. 6.3). 

4.3. Justification of the Changes for the Seismici Model 

The adjustments made increase the overall agreement between the Sun and the solar 
model on the c s profile (see Fig. 3). The improvement is less obvious when checking at the 
density profile but the progress is real too. To reach this improvement in the core and the 
radiative zone, we restricted the modifications to two parameters, S(0) p - P and Z . With 
such a restriction, it can be shown that we need to increase both these quantities to cancel 
(at least reduce) the discrepancy on c s . 

First, G&N93 work out (Z/X) s = 0.0245 ± 10% at the solar surface; this ratio has 
an impact on Z , whose value is larger than the present photospheric Z value because of 
the microscopic diffusion. Let us suppose we want to cancel the sound speed difference by 
increasing Z and reducing (or not changing) S(0) p - P : to keep a (Z/X) s ratio within its 
error bar, we must limit the raise of Z to ps 5.5%. This raise is lower than the one needed to 
cancel 5c 2 s /c 2 s . Thus, if we increase Z , we also need an increase of S(0) p - P to get Sc 2 /cj. = 0. 
Second, let us suppose we diminish Z : to cancel the sound speed difference in the core, we 
need to increase S(0) p - P beyond its upper error bar (equal to 2.0%). 

Thus, the error bars on S(0) p - P and on (Z/X) s imply an increase of both parameters 
to cancel the discrepancy on c s below 0.6R Q . 

4.4. Non Uniqueness of the Seismici Model: Alternative Seismic Models 

Of course, the seismici model is not unique: 
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First, we could have changed other nuclear cross sections instead of S(0) p - P . Yet, 
Turck-Chieze et al. (2001a) highlight the lack of sensitivity of the sound speed to the other 
nuclear reactions, like 3 He( 3 He, 2p) 4 He, 3 He( 4 He, 7) 7 Be or the CNO bi-cycle. For instance, 
an increase of 25% for the first two reactions mentioned changes 8c 2 s /c 2 s by only k, 0.1 — 0.2%, 
and thus does not improve significantly the agreement with the Sun. On the other side, the 
sound speed is sensitive enough to the p—p reaction rate. Therefore it seems more appropriate 
to only adjust S(0) p - p : we build the best model with the most simple assumptions. On Fig. 
3, you can see that the density is a little bit more sensitive to these reaction rates than the 
sound speed. 

Second, with the accuracy reached on the sound speed, this quantity is sensitive to 
many physical parameters whose uncertainties can be quite large. It is primarily sensitive 
to S(0)p- p , the opacities, the heavy element abundances, the microscopic diffusion process, 
but it can also undergo changes under "secondary" parameters such as the solar age or the 
solar radius. 



4.4.I. The Seismic2 Solar Model 

Actually, a second seismic model is constructed from Btz by modifying the opacities 
Kross, reducing Z by 3.5% and increasing S(0) p - P (we also calibrate this model to R 2 = 
6.95866 x 10 10 cm). Its main advantage is to reach a (Z/X) s ratio equal to the one proposed 
by G&N93, while the one of the seismici model approaches the upper limit of the error bar. 
By lowering Zq we reduce the mean molecular weight in the central part of the Sun, and the 
opacities as well. On the other side, we raise K ross following Brun et al. (1998), by simulating 
an increase of about 7.5% in the C, N, and O opacities. This raise could be attributed to 
uncertainties in the bound-bound and bound-free processes, i.e. to an intrinsic error in the 
opacity calculations. The change in the opacities of the CNO elements produces an increase 
of about 1.5% in n ross in the solar core, and about 4.5% below the convection zone. This 
raise compensates for the decrease due to the reduction of Z , so that we get roughly: no 
opacity change in the core and an increase of about 1.5% at O.6i? . We just operate this ad 
hoc adjustment below O.6i? to diminish the discrepancy on c s at the BCZ. As for the first 
seismic model, we also need an increase of S(0) p - P : by 1.3%. This seismic2 model produces 
a sound speed profile similar to the one of the seismici model, and even better at the BCZ 
(see Fig. 4). 

However, for the rest of the paper, we just consider the seismici model to derive the 
neutrino fluxes and to test the magnetic field: this model contains no ad hoc adjustment, less 
parameters were changed & the error bars on S'(0) p _ p and Z (through the (Z/X) s ratio) 
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are more or less well known while the uncertainties on K ross remain unknown. However, the 
great effort done by laser measurements and the theoretical successes of the Livermore group 
lead to specific uncertainties on K ross of the order of 5%. 

4 -4 -2. The Opacities and the Change in Metal Abundances 

Just a few words about these opacities: the OPAL95 tables we use to interpolate the 
K ross values were calculated assuming a G&N93 mixture. Along the solar evolution the metal 
composition in the core changes: because of the CNO poly-cycle and the microscopic diffu- 
sion, the relative number fractions of the metals are modified (see Fig. 2 for the composition 
at 4.6 Gyr). This should have an impact on the opacities: we have to correct them as soon 
as the solar core composition is "far" from the G&N93 mixture (this occurs at an age of, 
roughly, 100 million years). Starting from the seismici model, we compute a very similar 
model but with a correction for K ross below 0.15 R Q to take into account the change in the 
composition. Actually, this correction is very small: from —0.88% at 0R Q to —0.5% at 
0.15i? Q , for the Sun at 4.6 Gyr. It has a minute impact on the sound speed profile: it just 
increases b~c 2 s /c 2 by 0.0005 in the core! However, with the precision we reach on c s , we can 
test such a small effect. This latter leads to that, instead of increasing S(0) p - P by 1% for 
the seismici model, an increase by only 0.75% is enough. 



4-4-S. The Impact of the Age: Seismic^ 

Finally, we also check the sensitivity of c s to the solar age (£©). A recent work by 
Dziembowski et al. (1999) used helioseismology to determine t & . Depending on the method 
they use, they derive very different results. However, using the small frequency separation 
which they claim to be the most accurate measure, they conclude that t & = 4.66 ± 0.11 
Gyr (the lifetime on the PMS must be added). Starting from Btz, we compute the seismic3 
model: those are basically the same models except that t & = 4.735 Gyr for the second one, 
the tachocline parameters are changed (as for the previous seismic models), and we calibrate 
at i? 2 . The augmentation by only 2.9% of t & produces a rather important change in the c s 
profile and reduces the discrepancy with the Sun, eventhough the reduction is less significant 
than with the other seismic models (see Fig. 4). The seismic3 model shows that, supposing 
we under-estimate the solar age, an increase of S(0) p - P and Zq less important is needed to 
retrieve the flat Sc 2 /c 2 profile of the seismici model. Eventhough this does not challenge 
the need to raise both S(0) p - P and Z , this reduces the amplitude of the increase. On the 
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contrary, a solar age less than 4.6 Gyr would favor a larger increase of the two parameters. 
It is noticeable that the density profile greatly disfavors the seismic 3 model (see Fig. 4, right 
figure, panel b), unlike the sound speed profile. 

4.5. Constraints Suggested by the Seismic Models 

In the previous sections we emphasize the sensitivity of c s to a large set of parameters 
used to work out different seismic models: that makes it difficult to draw firm conclusions 
concerning the changes in just a few parameters, since these modifications cannot be proved 
to be the unique solution. However, under the assumption that the major uncertainties in the 
solar models (for the core and radiative zone) are due to some reaction rates, the metallicity 
and the opacities, our seismic models shed new lights on these physical quantities. 

First, the different seismic models we computed all favor a slight increase in the S(0) p - P 
factor: from 0.75% to 1.3% depending on the model. The strong influence of the p — p 
reaction was known for a long time (Turck-Chieze & Lopes 1993) and a raise was also 
recently proposed by Antia & Chitre (1998). However, the increase we work out is less than 
theirs (due to the data improvement) and could be even less if the solar age turned out to 
be slightly larger than 4.6 Gyr (on the contrary, a lower age implies a larger increase of 
S'(0)p_p). Incidentally, this result confirms that the calculations of S(0) p _ p are quite good 
despite their purely theoretical base. 

Concerning Z , an increase by 3.9% is strongly favored by the seismici model, while the 
seismic 2 model shows that an increase of K ross combined with a decrease of Z yields the 
same result in terms of c s . Since either solutions are acceptable, we cannot favor the increase 
or decrease of Z . Opacities and heavy element abundances are closely related, and part of 
the error bars on n ross are due to the uncertainty on Z. 

However, this change by almost 4% in Z is far from trivial if you consider that the 
microscopic diffusion changed Z by only 10% since the beginning of the hydrogen burning. 
Actually, the need to increase Z Q — in the seismici model — could be the manifestation of 
some forgotten hydrodynamic phenomenon. The use of both the density and rotation profiles 
in the core may help to solve this point. Gravity modes may also be helpful. For now, we 
only have a static view of the radiative region (contrary to the convective one). For instance, 
we cannot introduce the history of the angular momentum in the stellar equations: there- 
fore, it seems that only hydrodynamical simulations will be able to reproduce the real Sun. 
Nevertheless, the present analysis tends toward a reduced effect on the neutrino prediction 
when we do not correctly simulate the dynamical phenomena. Thus, this static view of the 
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Table 1. Main Features of the Solar Models 
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0.025 


0.025 


0.025 


0.025 


N (//Hz) 


25 


105 


45 


45 


45 


45 


fi (nHz) 


415 


430 


430 


430 


430 


430 


BCZ (i? ) 


0.7142 


0.7115 


0.7144 


0.7121 


0.7137 


0.7122 


71 Ga (SNU) 


127.1 b 


128.0 


125.9 


127.9 


135.1 


129.8 


37 C1 (SNU) 


7.04 b 


7.47 


7.11 


7.47 


8.88 


7.81 


8 B (10 6 cm- 2 s- 1 ) 


4.98 b 


4.98 


4.71 


4.98 


6.04 


5.23 


Atm c 


H 


K 


H 


H 


H 


H 



a This model was run with a time step larger than the one required to derive a realistic 7 Li 
abundance. Therefore, this quantity is non relevant here. 

b The neutrino fluxes for this model were derived with an old value for the capture rate of 
8 B neutrino by CI (1.06 x 10~ 42 cm 2 instead of 1.14 x 10" 42 cm 2 ), and for the S 17 (0) factor 
(19 eVb instead of 18.3 eVb). The impact of updating the ^(O) factor is to reduce the 8 B 
neutrino fluxes by about 3.5%. 

c The atmosphere model: H for Hopf, K for Kurucz5777 
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Sun remains an astrophysical problem that could have a larger impact on other stars. 

Concerning t Q , the sound speed and density profiles change in an opposite way: while 
a raise in t Q reduces the discrepancy on c s below 0.2R Q , it increases it for p. If we refer 
to the density profile, which is more sensitive, it seems that t =4.6 Gyr is a satisfactory 
value for the solar age (the value needed to best reduce the discrepancy on p below O.6i? is 
4.55 Gyr). Moreover, an augmentation of t is strongly disfavored by the density, while the 
sound speed does not rule it out (but does not strongly support it too). 

Finally, it seems that calibrating the solar models to a radius slightly smaller than 
the usual one (by 50-125 km, this value depends on the atmosphere model) improves the 
agreement between the Sun and the models above 0.9 R Q . 



5. THE OSCILLATION FREQUENCIES 

An important result of the solar models is the computation of the oscillation frequencies 
for the p and g modes. We use the seismici model: it was derived using the MHD EOS 
(Mihalas, Dappen & Hummer 1988) and the k5777 T(r) law (derived from a Kurucz's model) 
for the atmosphere model (with a reconnection at r = 20 following Morel et al. (1994)). This 
is an improvement in comparison with our previous Btz model using the T(r) Hopf's law, 
with the CEFF EOS (Christensen-Dalsgaard & Dappen 1992). 

The Figure 5 shows the difference — scaled by the usual Q n> i factors — between the 
theoretical and the observed p-mode frequencies, up to £ = 200, versus the inner turning 
point r t expressed as: 

w 2 + 4nGp _ L 2 
c s ' t 

w is the mode frequency, G the gravity constant and L = (I + 0.5) with £ the mode 
degree. We have taken the self-gravity of plane acoustic waves into account. 

As usual, there is a very good agreement for frequencies smaller than 2 mHz and a worse 
one at larger frequencies, mainly due to the difficulty to model the turbulent surface. The 
Table 2 gives the frequencies of the low degree g modes. 

The Figure 6 details the quality of the seismic data we use by showing the weighted 
difference between calculated and observed frequencies, as a function of the internal turning 
point. The surface effect has been removed by a polynomial fit of the general trend above 
2.2 mHz. We separate the modes in two ranges: in the first two figures we draw modes 
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with n > 14 and in the third one, those with n < 14. It appears clearly that the use of the 
higher frequency domain is not disfavored by the general trend previously observed. The 
real difficulty comes from the stochastic excitation and the mode correlated lifetime. This 
leads to a bad determination of the frequencies after 8 months of measurement (Fig. 6a). 
The situation is improved by a longer integration time (here 1290 days) (Fig. 6b) but it is 
obvious that the comparison computed/observed frequencies is better when we access to the 
low frequency range (n < 15 or v < 2.2 mHz). In this case, the modes have a longer lifetime 
and an insight into the very central core (which is important for high energy neutrino fluxes) 
could be obtained with a large increase of the sensitivity. Therefore, the inversion of the 
sound speed is better determined (even if the use of the sole low order modes reduce the 
radial accuracy), as for the density profile. This is natural since we also reach the modes 
that have a mixed gravity and acoustic character. 



6. SEISMIC PREDICTED NEUTRINO FLUXES 

6.1. The Solar Neutrino Predictions 

Since the neutrino production is not pointlike, the precise calculation of the expected 
emitted v e flux requires to know where the neutrinos are produced. The Figure 7 recalls the 
production zones for the p-p, pep, 8 B, 7 Be, 13 N, 15 and 17 F neutrino fluxes. As it clearly 
appears on the figure, most of them depend on the very central part. Now we recall that 
the first radius we have for the sound speed is at 0.07 ± 3.6%i? (Table 1 of Turck-Chieze 
et al. 2001b), and only gravity modes may improve this situation. 

Neutrino predicted fluxes have been calculated for the three seismic models (see Table 1). 
We may note that they do not differ significantly. Moreover, they are only slightly different 
from our previous predictions in Brun et al. (1999) or from the recent solar models (e.g. 
Bahcall et al. 2001). This is logical since the last modifications we introduced are minor. In 
Turck-Chieze (2001), we recall the progress made to stabilize these fluxes. Detailed neutrino 
capture predictions are shown for the seismici model (see Table 3). 

The present predictions are of great interest (in comparison with those obtained ten years 
ago) since they include helioseismic data that validate the updated physics at the needed 
level (see Table 2 of Turck-Chieze et al. 2001a). The emitted neutrino fluxes we predict are 
not only theoretical but also partly deduced from precise seismic "observations" of the solar 
core, and they confirm the discrepancy between predictions and detections, favouring the 
neutrino flavor transformation. The next step is to sum the fluxes of these different flavors 
received on Earth. Such data are now available thanks to the SNO experiment. 
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Table 2. G-mode Frequencies with the Seismic i Model 
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95.47 
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-13 


136.49 
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-15 
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170.87 
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126.96 


4 


-9 


179.62 


5 


-10 


192.39 


6 


-12 
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6.2. Comparison of the 8 B Emitted Flux with the SNO & Superkamiokande 

Results 

Ahmad et al. (2001) announced the first results of the SNO collaboration on charged 
current (CC) and elastic scattering off electrons (ES) reactions. 

The CC reaction is only sensitive to electron-type neutrinos, while the ES reaction, also 
measured in the Superkamiokande experiment, is sensitive to all flavors. The measured 8 B 
neutrino fluxes for the CC and ES reactions are: 

$g^ (z/ e ) = 1.75±0.07(stal)+°;^( sys .)l0 6 cm- 2 s" 1 

$sivoK) = 2.39 ± 0.34(stat.)+°;^(sys.)10 6 cm" 2 s" 1 

The present statistical error on the elastic scattering is still large and the sensitivity to 
the v^, v T or anti v T (only 16% in comparison to 100% for u e ) is low. Therefore, it is not 
possible yet to properly extract the different flavors from the SNO experiment alone, despite 
the difference in the two measured fluxes. Fortunately, the ES result is in agreement with 
the measure from the SK experiment (Fukuda et al. 2001): 

$f£(i/ x ) = 2.32 ± 0.03(stat.)l{j;g?(sys.)10 6 cnr 2 s- 1 

This flux is equal to ~ 47 — 48% of the seismici model 8 B neutrino prediction (Turck- 
Chieze et al. 2001b): 

&{y x ) = 4.98 ± 0.73 x lO^nrV 1 

It is widely known that the 8 B neutrino flux is the most difficult to predict and is very 
dependent on the physics of the core. A fundamental improvement achieved by the seismic 
model is the reduction of the error bars in the neutrino prediction, through the reduction of 
the uncertainty on the p — p reaction rate, the (Z/X) s ratio... This point has already been 
discussed in the Table 2 of Turck-Chieze et al. (2001b) that gives the detailed uncertainties 
on neutrino predictions. We point out that the neutrino fluxes derived here are not only 
"observational" (to a certain extent) but are also affected by reduced uncertainties. They 
also include the recent rejection of several astrophysical solutions to the neutrino puzzle 
proposed in the past. 

Nevertheless, the 5(0)17 factor has no impact on the structure and was badly determined 
more than 4 years ago (see Turck-Chieze 2001). For instance, Adelberger et al. (1998) propose 
5(0)17 = I9I2 e V b. But the experimental measurements have been largely improved. In 
our models, we use the value of Hammache et al. (1998) of 18.5 ± 1 eV b. Unfortunately, the 
recent result of Junghans et al. (2002) of 22.3 ± 0.7 eV b, is only marginally in agreement 
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with the previous one. Thus the uncertainty on the present predicted neutrino flux, including 
the one of the seismici model, could be a bit larger today. 

Actually, the major result of the last year is the one obtained by the SNO collabora- 
tion when they add the estimated number of detected muon and tau type neutrinos to the 
electronic ones (using the SK results); they find: 

$(z/ x ) = 5.44 ± 0.99 x 10 6 cm- 2 s _1 

a result very consistent with the previous value validated by the present helioseismic data. 

This strongly favors or even "proves" the existence of neutrino oscillations: a part of 
the v e must be converted into and u T , as far as new seismic results do not lead to other 
conclusion. Actually the present study stays in the framework of a static solar core, with 
classical phenomena. This representation is compatible with the present seismic results but 
we still need to extract the rotation and magnetic field in the radiative zone to check this 
assumption. 

6.3. The Neutrino Oscillation Related Quantities 

6.3.1. The Electron Number Density 

Following this framework, a part of the neutrino oscillations could be explained by the 
well-known Mikheyev-Smirnov-Wolfenstein effect (MSW effect, see e.g. Mikheyev & Smirnov 
1986): an electron type neutrino may undergo a resonant oscillation in the Sun and then 
be converted into a muon or tau type neutrino. This effect assumes that the neutrinos have 
masses and that the flavor eigenstates are different from the mass eigenstates. 

For the sake of simplicity, we suppose two flavors [y e and v^) and two mass eigenstates 
[y\ and v-i). We assume that v\ is close to v e , and u 2 is close to v^. When an electron 
type neutrino is created in the solar core, the flavor eigenstates are almost "propagation 
eigenstates", because of the high electron number density (n e ). But during the propaga- 
tion from the neutrino toward the solar surface, n e decreases and the flavor states are no 
more propagation eigenstates: the neutrino state vector starts to oscillate around the new 
propagation eigenstate. This latter changes with the further decrease of n e , and gets closer 
to u 2 (initially, it was close to u e ). Depending on the variation of n e , this change is more 
or less rapid and said adiabatic or not. In the case of an adiabatic change, the probability 
for an emitted electron type neutrino to be detected as a muon type neutrino at the solar 
surface is rather simple to calculate. In the non-adiabatic case, the conversion probability is 
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more tricky to work out. In both cases, the electron number density is needed with a high 
accuracy all along the solar radius, to compute the conversion probabilities. 

Therefore, we derive n e for the seismici model. To do so, we use the precise mass 
fractions returned by CESAM: we know the mass fractions for the 1 H, 2 H, 3 He, 4 He, 7 Li, 
7 Be, 9 Be, 12 C, 13 C, 14 N, 15 N, 16 0, & 17 chemical elements, plus an extra element (with 
A = 28 and Z = 13), as a function of the fractional radius. We just have to derive the 
number fractions of all these elements, and multiply these number fractions by the electron 
number of the corresponding chemical element. By adding all the terms we obtained this 
way, we determine the electon number density as a function of the radius. The result is 
drawn in the upper part of Fig. 8. 

6.3.2. The Neutron Number Density 

If the neutrino has a magnetic moment (either a dipole or/and transition moments), it 
might interact with the solar magnetic field. Provided that the magnetic moment and the 
magnetic field are large enough, this field could flip the spin of the neutrino: a left-handed 
neutrino could become right-handed. Moreover, the possible flavor transition magnetic mo- 
ments could result in a spin-flavor precession: the neutrino could change both its chirality 
and its flavor. This double precession could be matter-enhanced through the interactions of 
neutrino with the electrons, protons and neutrons: this is the Resonant Spin-Flavor Preces- 
sion process (RSFP, e.g. Lim & Marciano 1988). To compute the conversion probabilities for 
the RSFP, the neutron number density n n is needed (see the lower part of Fig. 8): we derive 
it the same way as the electon number density, but instead of using the electron number of 
each chemical element, we use its neutron number. It is important to note that both the 
MSW and RSFP processes could cohabit inside the Sun. 

The detailed n e and n n profiles are available with the seismici model on the web site 
whose address was previously mentioned, in order to calculate quantities related to the 
neutrino oscillations. 



7. UNSOLVED PROBLEMS IN THE SEISMIC MODELS 

Despite the overall agreement in the sound speed between the Sun and the seismic 
models below 0.6i? Q , two regions of our star remain poorly described: the tachocline region 
and the upper layers, both are expected to undergo dynamic effects. Moreover the central 
rotation law is not taken into account in this analysis (Garcia, R. A., Couvidat, S., Turck- 
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Chieze, S. et al., in preparation) . 

In the superadiabatic region it is well known that the turbulent pressure becomes non 
negligible compared to the gas pressure. Moreover, the tachocline and the upper solar layers 
are shear layers in which rotation rates change rapidly (e.g. see the rotation curve of Howe et 
al. 2000). A ID stellar evolution code cannot afford an efficient treatment of these dynamic 
regions. Neither the rotation of the Sun nor its magnetic field are taken into account, 
whereas it is widely thought that the tachocline is the base of the magnetic dynamo process. 
Of course, the neutrino emission and the solar core physics are rather insensitive to the 
tachocline and beyond, but the neutrino behaviour may depend on these layers. Progresses 
are needed to connect interior phenomena related to these neutrinos, to what is observed 
out of the star. If the neutrinos have a magnetic moment, they could interact with the solar 
interior magnetic field. In the next section we focus on the solar large-scale magnetic field: 
we add magnetic pressure and derive new solar models. 



7.1. The Magnetic Field Tested by the Sound Speed 

The goal of this analysis is to test the sensitivity of c s (r) to the magnetic field and 
constrain the amplitude of this field. We also check (and use) the sensitivity of the density. 

Basically, the presence of a magnetic field changes the wave velocity in two ways: 

• it changes the gas pressure because of the hydrostatic equilibrium; 

• it adds a part of the Alfven speed v a = B 2 /(Airp) (in cgs units) to the wave velocity. 

Thus, any magnetic field should be imprinted in the "magneto-acoustic" wave velocity. 
Of course, the 3D structure of the field is of great importance since the angle between the 
field lines and the seismic waves determines the way these latter are accelerated: the wave 
velocity is no more an isotropic quantity. Unfortunately, we cannot account for the field 
structure with a ID stellar code, and so we just add a magnetic pressure term P mag in the 
stellar structure equations: P mag = B 2 /(8tc) (in cgs units). Therefore, the basic hydrostatic 
equilibrium equation is changed into: 



dPgas _ GM dP ma g . . 

dm Anr 4 dm 

The main problem is to choose an appropriate magnetic field B(r) for the solar interior, 
since little is known about the inner field. Roughly speaking, it seems that large scale 
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magnetism may be important in three regions: the radiative zone, the tachocline, and just 
below the solar surface. 



7.2. Origin of a Large Scale Magnetic Field 

Parker presented the hydromagnetic dynamo theory for the Sun in 1955. The a — uj 
kinematic dynamo needs both differential latitudinal rotation (parameter a;) and turbulent - 
cyclonic — movements (parameter a). The differential rotation produces a toroidal magnetic 
field from a poloidal one, while the cyclonic movements slow down the lift and twist of 
this toroidal field which has enough time to get amplified. The lift of the field is due to the 
magnetic buoyancy, while the twist results from the Coriolis force. Both these actions induce 
a poloidal field from the toroidal one, thus restoring the initial field. This basic picture has 
been widely reviewed as different dynamo processes were developped and improved. For 
instance, let us mention the Babcock-Leighton dynamo models that regenerate the poloidal 
field by the eruption of the toroidal field at the solar surface (e.g. see Durney 1997 for a 
modified Babcock-Leighton dynamo analysis). It is now commonly believed that the seat 
of the (main) solar dynamo process is the tachocline. A small seed poloidal field — a legacy 
of the PMS evolution — in the radiative interior should have given birth to a magnetic 
field with poloidal and toroidal components. Although this scenario seems widely accepted, 
the situation is rather confused when going into details. The seat of the poloidal field 
regeneration is not located, it could be the tachocline (as is induced by the "Parker" dynamo 
process), or just below the solar surface (as proposed by the Babcock-Leighton dynamo). 
Moreover, it is unclear whether or not the toroidal part of the field prevails in the solar 
interior: the surface activity proves that the toroidal field is larger than the poloidal one in 
the upper layers of the Sun, but only a few clues of what happens in the deep interior are 
known. We just consider toroidal fields in this paper. 

7.3. Simulated Magnetic Profiles 

Following Dzitko (1995) we simulate fields as: 

B<P = a ( r )~^Q P k( C0Sd ) e <t> ( 3 ) 

With the spherical coordinates (r, 9, 0). Pk(cos9) is a Legendre polynomial of degree k. 
For argument's sake, we assume k = 2, meaning that the field is quadrupolar (in ac- 
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cordance with the surface-magnetism manifestations). The presence of a magnetic pressure 
modifies in a direct way the wave velocity: we must take into account the Alfven speed. 
Since we are only interested in the radial velocity and a toroidal field is perpendicular to 
the radial direction, the new wave velocity is \J & s + v\ (on the figures, we continue to note 
this quantity "c s " for convenience, eventhough it is no more the sound speed, strictly speak- 
ing). For the function a(r), two profiles are considered, following Gough & Thompson (1990). 

First, to simulate a magnetic field in the radiative zone, we choose: 



With K x = (1 + A)(l + 1/\) X B , r = O.712i? being approximately the BCZ and 
A = 10r + 1. Bq — the highest intensity of the field — is set to different values (see Table 
4). The addition of the magnetic pressure is made when the Sun enters the ZAMS. The 
intensities of these fields are very high, compared to the prediction of Gough & Maclntyre 
(1998) who claim for a (poloidal) field with an amplitude of ~ 0.1 T in the deep interior. 
Moreover, the virial theorem rules out fields larger than w 10 4 T (Mestel & Weiss 1987). 
The maximum values of the P mag /P g as ratio (hereafter f3~ l ) are available in Table 4, for all 
the models discussed here. We could also give the v a /c s ratio, but (3 and this ratio are closely 
related: v a /c s = a/2/ {fi f\). Both contributions of the magnetic field to the change in c s - 
the indirect change through the modification of the solar structure, and the direct change 
through the addition of v a — depend on the (5 values. 

Second, we also utilize as a(r) profile: 



With d the half-width of the zone with a magnetic field, and tq its center. To simulate 
a magnetic field in the tachocline, parameters are set to: d = O.O2i? and the radius r 
of transition between radiative and convective zones varies along the solar evolution. We 
set B to 30 T according to Antia et al. (2000) for the seismiciB 2 model, and 50 T for the 
seismiciB 2 i model. 

Finally, this profile is also used to simulate a possible field in the upper solar layers. 
Such a field was hinted by Antia et al. (2000) from an analysis of the Global Oscillation 
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Table 3. Neutrino Fluxes of the Seismici Model 
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Capture 


Predictions for 
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71 Ga: 


128.1 ±8.9 SNU, 


37 C1 : 7.48 ± 0.97 SNU, water experiments 


: 4.98 ±0.73 


x 10 6 cm^s" 1 





a ' b flux emitted by the Sun, and received on Earth, in units of neutrinos x cm~ 2 s _1 

c in SNU, Solar Neutrino Unit: 1 SNU = 10 -36 captures per second and per target atom 

d The error bars on the predictions are derived from table 2 of Turck-Chieze ct al. (2001b) 



Table 4. Models with Magnetic Field 
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a radius at which P mnn is maximum 
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Network Group (GONG) and MDI data. To test a field anchored at 0.96_R Q , parameters 
are set to r = O.96i? and d = O.O35i? . B is set to 2 (the value proposed by Antia et al. 
2000) and 3 T. The fields at the BCZ and in the upper layers are both added at 85 million 
years, when the solar core is no more convective. 

Since all the physical quantities considered here are radial quantities, we average the 
magnetic pressure P ma g(r,9) over 9. 

7.4. Impact of the Magnetic Pressure on the Solar Model: c s , p, and the 

Neutrino Fluxes 

The three different magnetic pressure profiles are drawn on Figure 9. The addition 
of magnetic pressure in the radiative zone following the seismiciBi model induces a great 
change in the thermodynamic quantities (P,T...), especially in the c s and p profiles (see Fig. 
10 & 11). On the contrary, the impact of a field 10 times smaller following the seismiciB 13 
model is minuscule (see Fig. 10 and 11). The precision we have on the solar sound speed and 
density rules out a magnetic field with such a profile and an intensity as large as B$ = 10 4 
T. Actually, we can put an upper limit for a (toroidal) magnetic field in the radiative zone 
of about 3 x 10 3 T. As soon as B < 10 3 T, c s and p are not sensitive enough to P mag and 
we cannot draw any conclusion about the likelihood of a field like the one of the seismiciB 13 
model. We show in Table 1 the impact on the neutrino production of the seismiciBn model. 
The fluxes are slightly larger than the ones deduced from the adjustment in the physics 
of the seismic models. The result obtained for the 8 B flux: 5.23 x 10 6 cm _2 s _1 remains in 
agreement with the present SNO result. 

We face the same problem when adding P mag at the BCZ and in the upper layers (see 
Fig. 10). It is even more difficult to draw any conclusions: the sound speed and density are 
not sensitive enough to such modifications. With the current accuracy we have on c s and p, 
it is only possible to state that a field in the tachocline can reach an amplitude as large as 
50 T without perturbing the sound speed profile. We conclude the same for a toroidal field 
anchored at O.96-R and as large as 3 T. Given the minuscule impact on c s and p, we did 
not draw the 5c 2 s /c 2 s and dp/ p profiles obtained with the fields at the BCZ and in the upper 
layers. 

This section on the magnetic field confirms that c s and p are only sensitive to the f3 
ratio. With the different models we computed, each one with a different magnetic field profile 
and/or intensity, we can conclude that only the large scale fields with larger than, at 
least, ~ 3 x 10~ 4 , impact on the sound speed and density profiles. The firm conclusion to 
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draw is that an intensity as large as 3 x 10 3 T can be ruled out for a in-radiative- zone toroidal 
field. Despite the great accuracy we reached on this quantity the sound speed is not suited 
to the determination of the large scale magnetic features of the Sun: many physical processes 
c s is quite sensitive to, are still affected by large error bars, and a potential magnetic field 
looks like "background noise" compared to these processes. The same conclusion is valid 
for the density profile. Yet the use of these two quantities might be more promising for 
constraining an upper-layer field, since the "weakness" of P gas near the solar surface makes 
c s more reactive to weaker fields. Unfortunately, this part of the Sun is also badly modeled 
due to the lack of turbulence in the present stellar equations. 

Concerning the neutrino puzzle, none of the reasonable models including a central mag- 
netic field (as models seismiciBn, seismiciB 12 , & seismiciB 13 ) greatly modify the neutrino 
flux predictions, except the ruled out seismiciBi model (it increases the 8 B neutrino flux by 
more than 20%). With the upper bounds we have for the large scale magnetic field, it seems 
that this field has no impact on the neutrino emission, or only a very slight one. 

7.5. About the Relation Between Magnetic Field and the Neutrino Transport 

So far, we just considered mean static fields. We concluded that these fields probably 
impact only slightly on the neutrino production. However, the solar magnetic field is expected 
to be highly variable, and may locally be of large intensity: in the convective zone, the field 
is expected to be concentrated in flux tubes. Actually, the very idea of a mean field seems 
to be fallacious, at least in the external part. A variable and local magnetic field could act 
on the neutrino transport, since this field is expected to be by far larger than an hypothetic 
mean field. 

We know from our previous analyses that a field with P ma g/Pgas = 10~ 4 changes the 
gas pressure by at most 5P/P = —1.65 x 10~ 4 . It changes the sound speed by Sc^/c^ = 
—0.74 x 10~ 4 , and the density by Sp/p ~ — 1 x 10 -4 (the minus sign means that when you add 
a magnetic field, you reduce the density and the sound speed). Of course, a change in the 
density means a change in the electron number and neutron number densities. Depending 
on the flux tube magnetic profile, and on the intensity of this field, such changes can impact 
on the neutrino transport: the appearance of a flux tube on the neutrino path may flip the 
neutrino spin, the sudden change in the electron and neutron number densities may act on 
the neutrino oscillations... 

To correctly account for all these phenomena, we will have to use the results of magneto- 
hydro-dynamic simulations of the solar magnetic tubes. Thus we will get access to the way 
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the electron and neutron number densities vary inside the tube, and the magnetic profile of 
this one. With a ID stellar code, we cannot reproduce the complex magnetic structure of 
the Sun, but we can begin to work out some estimates, very near the surface where the local 
magnetic field is known and dominates the gas pressure. We are ready to begin this study 
to estimate the impact on neutrino oscillations. 

8. CONCLUSION 

Thanks to the most recent seismic data from the GOLF and MDI instruments aboard 
SoHO, precise sound speed and density profiles were derived down to 0.07 R , that enabled 
us to build solar seismic models. The main goal was to obtain models in total agreement 
on the sound speed with the real Sun in the regions where no dynamic effects have been 
observed yet. The two first seismic models proposed perfectly fulfill this requirement below 
0.6 R Q : it seems that further improvements in solar modelling require 3D codes to correctly 
account for the dynamic processes in the Sun. 

With the seismic models, we derive frequencies for acoustic and gravity modes and neu- 
trino emitted fluxes: these fluxes are no more the results of purely theoretical considerations, 
but they also include seismic data. They are known with a better precision. Recent releases 
from the SNO collaboration announce a 8 B neutrino flux received on Earth very close to 
the one predicted here. Of course, this is a great breakthrough in the understanding of 
the neutrino puzzle. The seismic models considered so far are the result of long years of 
improvements in both the seismic data and the solar physics. The sound speed is now quite 
sensitive to the p — p reaction rate, the metallicity, the opacities, the solar age and even to 
some magnetic fields. This is a very powerful tool that might be useful to constrain the large 
scale solar magnetic field, provided that some further improvements are realized on the other 
physical parameters that have a large impact on c s . It seems that the density profile should 
be closely monitored since it is more reactive than c s to the changes in several physical pro- 
cesses. The present density profile recently derived needs to be confirmed by an enhanced 
number of low frequency modes. An important feature of the seismici model is to favor a 
slight increase of the p — p reaction rate by ~ 1% and the initial metallicity by about 3.5%. 
Actually, it is rather difficult to choose between an increase of the metallicity and a change 
of both the opacities and the metallicity. It also seems that the models must be calibrated at 
a radius smaller than the usual one (by 50 — 125 km, depending on the atmosphere model), 
to improve the agreement with the Sun in the upper layers. This proposition needs to be 
cautiously considered, since a complex physics is present in the external layers, including 
turbulence and magnetic field. Finally, the solar age 4.6 Gyr is a good compromise between 
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an older Sun favored by the sound speed profile and a younger one favored by the density 
profile. However, the density strongly disfavors an increase of t Q . 

We put constraints on the magnetic field in the solar core and show that the fields 
simulated here do not deteriorate the agreement with the SNO observations, supposing 
three neutrino flavors. We compute the different quantities useful to deduce the neutrino 
oscillation parameters, in the framework of classical models. We also begin to estimate to 
what extent these classical models are representative of the real Sun. Many progresses are 
still needed in the solar models — more specifically for the upper layers — but the high quality 
of the seismic data combined with the improvement in the physics make us already achieve a 
good result in solar modelling: the seismicx model seems very close to the real Sun, and the 
seismiciBn model is also interesting to consider due to the presence of magnetic pressure. 
Contrarily to the usual seismic models, ours were obtained with a classical stellar evolution 
code. Nevertheless, as previously mentioned, as far as the very internal rotation profile is not 
included in such a study, new surprises may appear that invalidate the classical approach as 
it does for the solar convective zone. 
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Fig. 1. — Sound speed and density inversion using a standard Saclay model (Brun et al. 
1999). In these plots, panels a) and b) are obtained from the GOLF modes of 1=0-2 above 1 
mHz (Bertello et al. 2000a, b) and MDI medium-l data (1=3-250); the total number of modes 
is 2183. Panels c) and d) are obtained by adding three modes provided by Bertello et al. 
(2000b) (1=0, n=3,5,6) to the previous combination of p modes. 
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Fig. 2. — Upper figure: the contribution of the heavy elements to the opacity (plain curve). 
Superimposed is the profile of the Rosseland opacities n ross (dot-dashed curve). Lower figure: 
number fraction of three elements related to the metal number fraction: the carbon (plain 
curve), the nitrogen (dashed curve), and the oxygen (dot-dashed curve). Solar composition 
at 4.6 Gyr. 
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seismic^ model (right panels) (see text). 
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Fig. 6. — Difference between the p mode frequencies of the seismici model and the observed 
frequencies from GOLF, as a function of the internal turning point. In the upper figure: we 
have used the frequencies of Lazrek et al. (1997). In the middle and lower figures: Garcia 
et al. (2001) and Bertello et al. (2000a). For the upper and middle figures we show the 
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Fig. 7. — Neutrino production as a function of the fractional radius. In the upper figure, we 
have drawn the p-p (plain curve), 8 B (dot-dashed curve), 7 Be (dashed curve), and the pep 
(dot- dot- dot- dashed curve) neutrinos. In the lower figure, the 13 N (plain curve), 15 (dashed 
curve), and 17 F (dot-dashed curve) neutrino production is shown. For each neutrino "type", 
we have drawn (1/ F t ) (dF/dr) where F is the flux in s" 1 , r the fractional radius, and F t the 
total flux for this neutrino type (integrated over the entire Sun). 
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Fig. 9. — Magnetic fields simulated: in the radiative interior (plain curve), in the tachocline 
(dashed curve), in the upper layers (dot-dashed curve). The amplitudes of the fields have 
been normalized. 
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Fig. 10. — Upper figure: difference in the gas pressure between the seismic\B\ and the 
seismici models (thin plain curve). Same figure for the seismic\B\\ (dashed curve), the 
seismic\Byi (dot-dashed curve), and the seismic\B\2 (dotted curve) models. Superimposed 
is the difference in the temperature between the seismiciB 1 and the seismic\ models (thick 
plain curve). Lower figure: Same figure for the seismic\B 2 model (thick plain curve), the 
seismiciB 2 i (thin plain curve), the seismic\B 3 (thick dot-dashed curve), and the seismiciB 3 i 
(thin dot-dashed curve) models. 
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Fig. 11. — Upper figure: difference in the square of the sound speed between the Sun and the 
seismic\B\ model (thin plain curve). Same figure for the seismic\B\\ model (dashed curve), 
the seismic\B\2 model (dot-dashed curve), and the seismic\Bi 3 model (dotted curve). The 
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